A;-essence in the DGP brane-world cosmology 
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We analyse a DGP brane filled with a fc-essence field and assume the fc-field evolving linearly with 
the cosmic time of the brane. We then solve analytically the Friedmann equation and deduce the 
different behaviour of the brane at the low and the high energy regimes. The asymptotic behaviour 
can be quite different involving accelerating branes, big bangs, big crunches, big rips or quiescent 
singularities. The latter correspond to a type of sudden singularity. 



I. INTRODUCTION 

One of the most puzzling discovery of the last years in 
physics is the current acceleration of the universe [1| . De- 
spite the huge efforts made so far to find a well motivated 
theoretical framework for this behaviour no consensus 
has been reached, even though a fine tuned cosmological 
constant is the simplest option to match the current ob- 
servations [l|-[j] . There are two main streams of thought 
that try to explain the late-time speed up of the universe: 
(i) modified theories of gravity on large scale, which by 
weakening the gravitational interaction on those scales 
allow inflationary universes (cf. Refs. 0t3) (ii) a dark 
energy component, corresponding to a new component on 
the cosmic pie of the universe, that violates the strong 
energy condition and therefore allows accelerating uni- 
verses (lOj . 

One approach to build a modified theory of gravity re- 
lies on extra-dimensions. The idea of extra-dimensions 
is quite old in physics and it dates back to Kaluza and 
Klein. On the last decades it has been invoked by string 
theory as an approach to unify the different interactions 
in nature and in particular as a new road to obtain a 
consistent theory of quantum gravity. On this modern 
approach not only the extra-dimensions plays a crucial 
role but also the branes [llj . In particular, within the 
context of brane-world models, our universe corresponds 
to a brane; i.e. a 4-dimensional (4d) hypersuface, embed- 
ded in a higher dimensional space-time dubbed the bulk. 
The simplest of these models are Randall-Sundrum mod- 
els corresponding to an ultra-violet modification of gen- 
eral relativity (GR) [12[ or the Dvali-Gabadadzc-Porrati 
(DGP) model corresponding to an infra-red modification 
of GR [13| . The latter is a promising approach, despite 
its shortcomings [14[ , to describe the current inflationary 
epoch of the universe. 

The DGP model contains two set of solutions, usually 
referred to as the self-accelerating branch and the normal 
branch. While the self-accelerating branch as it name 



indicates is accelerating in the absence of any kind of 
dark energy, the normal branch requires some sort of stuff 
to describe any inflationary era of the universe. On the 
other hand, the normal branch is free from the ghost issue 
present on the self-accelerating DGP solution. However, 
it is important to highlight that there is a duality between 
both branches. Indeed, we can obtain the homogeneous 
and isotropic solutions and investigate them in anyone of 
the two branches (see section HI]). In the present paper, 
we show how it is possible to get a set of accelerating 
branes of the normal branch by means of a fc-essence 
field embedded on the brane [la . 1 1GJ ] . A complementary 
motivation for our analysis is to classify the different kind 
of singularities, in particular those related to dark energy, 
that may show up on the brane. 

The paper is outlined as follows. In section II, we 
present the model we will analyse, i.e. a homogeneous 
and isotropic DGP brane filled with a fc-cssence field that 
evolves linearly with the cosmic time of the brane. At the 
end of this section we introduce the form invariant trans- 
formation, that preserve the form of the Friedmann and 
conservation equations, to show that there is a duality 
between the solutions of the normal branch and the self- 
accelerating one in the DGP model. This duality relates 
expanding and contracting solutions among themselves. 
In section III, we present the analytical solutions of the 
model for the normal branch and analyse the different be- 
haviour of the brane depending on the equation of state 
of the fc-field and its energy density. On the last section, 
we summarise our main results and conclude. 



II. THE DGP MODEL WITH A fc-ESSENCE 
FIELD 

We will analyse the evolution of a DGP brane [13[ filled 
with a fc-cssence field [15| confined on the brane. The 
energy momentum tensor reads [15| 

dF , s 

(2.1) 
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Here F = F(x) is an arbitrary function of the kinetic 
energy of the fc-field x = g** v V ^V v <f> and V{4>) is a po- 
tential. We will assume a spatially flat, homogeneous 



and isotropic brane, therefore, the modified Fricdmann 
equation reads [13] 



H 1 



-H 



(2.2) 



where r c is the cross-over scale and is related to the 
ratio between the 4-dimensional effective gravitational 
constant, k|, and the 5-dimcnsional gravitational con- 
stant of the bulk. Finally, e = ±1 stands for the two 
branches of the DGP model, e = — 1 corresponds to the 
self- accelerating DGP branch and e = +1 corresponds to 
the normal DGP branch [17J. For latter convenience we 
rewrite the Friedmann equation as 



ff = ST 



■1± 4 /1 



where 
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The energy momentum tensor is conserved on the brane 
and therefore 



p + 3H(p + p)=0. 



(2.5) 



In addition, it can be shown that the Raychaudurui equa- 
tion can be written as 
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From Eq. (|2.1|) . we obtain the energy density p and 
pressure p of the fc-field 



V(4>){F-2xF x ), p = -V{4>)F, 



(2.7) 



with equation of state 7 = (p + p)j p = —2xF x /(F — 
2xF x ). The evolution equation of the fc-field can be de- 
duced by substituting Eqs. ([2~77]l into Eq. (j2~5j) 

(F x + 2xF xx )4> + 3HFJ + ^{F - 2xF x ) = 0, (2.8) 

with ' = d/ d<p. 

For simplicity, we consider the fc-field evolving linearly 
with time [161 [l8j; i.e. <f> = <f>ot where <fio = constant. 
Then, «o = F — 2xF x | a! __ $2 is constant and Eq. (|2.8[) 
implies 
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The previous equation can be integrated with result 
V = Voa~ 3jo . We obtain the same potential deduced 
in [la] because the result is independent of the modified 
Friedmann equation of the brane. Similarly, we can show 
that p = po a_37o i where po = ctoVo- 



In order to analyse the dynamics of the brane, it is 
useful to introduce the following variable 



z= Pm_ 



(2.10) 



Therefore, the Hubble parameter can be expressed as 



H 



(2.11) 



1 Z 

370 z 

Notice that by specifying the variable Z, we can deter- 
mine the fc-essence field potential as V ex 1/Z. This 
last relation can be deduced by combining Eqs. (|2.9p and 
(|2~TTj) and <p = (f> t. Finally, by substituting Eq. ([2TT) 
in the modified Fricdmann equation (|2.3[) . we obtain 
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(2.12) 



There are two aspects of the normal DGP branch that 
deserves to be commented: (i) it requires some sort of 
dark energy, for example in the form of a fc-essence field, 
to describe the late-time behaviour of the universe, while 
the self-accelerating branch does not require any cosmic 
fuel to mimic the current acceleration of the universe 
B LUJ; (h) the normal branch is free of ghosts, a charac- 
teristic presents on the self-accelerating branch [14| . 

At this point, it is interesting to note that the two 
branches of the DGP model can be related by a form 
invariant transformation that preserve the form of the 
Friedmann and conservation equations. In fact, the 
change 



H 



-H, 



-e, 



P- 



- P 



(2.13) 
transforms the two branches of Eq. (|2.2[) between them 
and preserves the form of the conservation equation (|2.5p . 
This internal symmetry gives rise to a duality: a — > 1/a, 
which after integration implies H — » —H, between ex- 
panding and contracting universes [IJI- More precisely, 
if we know a solution a of the branch e corresponding to 
a fluid with energy density p and pressure p, that satis- 
fies the weak energy condition, then 1/a is a solution of 
the branch — e with the same energy density and pres- 
sure —2p — p, that violates the weak energy condition. It 
means that there is a symmetry between the normal and 
the self-accelerating branches. 

Along this paper, we will use the normal branch and 
the conclusion for the self-accelerating branch will be 
obtained by applying the form invariant transformation 
(l2~T3l) . 



III. THE NORMAL DGP BRANCH 

In what follows, we concentrate our analysis on the 
solutions of a brane filled with the fc-field introduced on 
the previous section for e = 1, corresponding to the nor- 
mal DGP branch geometry. There are four cases to be 
analysed depending on the sign of the equation of state 
70 and the energy density of the fc-field p. 



A. Positive energy 



For < p; i.e. < Z, the Eq. ([2~T2]) can be integrated 
analytically [20 ] 



Z± 
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(3.1) 



where i^ and ti arc constants. The solutions of Eq. 
(I2.12p for + and — will be written as Z + , Z~ respectively. 

If the fc-ficld mimics a cosmological constant; i.e. 70 = 
0, then Eq. (|3.1|) is easily satisfied as the lhs and rhs 
of Eq. (|3.ip vanish (notice that the energy density and, 
therefore, the variable Z ± are constants on this particular 
case). 

At high energy p m < p, i.e. Z ± < 1 (see Eq. (|2~TU1) ). 
Z can be approximated as 






(3.2) 



where we have made a rescaling of the cosmic time. 
Therefore, V ± ex (</> /</>) 2 . We kept the square root on 
the previous equation because it is crucial to fix the sign 
of the cosmic time of the brane for a given sign choice of 
70. At low energy p <C p m \ i- c - 1 "C Z^, it can be proved 
that 

37o 



Z^ 



ln(Z-) 
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(3.3) 
(3.4) 



Consequently, the fc-field potential on this regime ful- 
fils: V + ex <J}q/4> and V~ oc cxp[37o</>/(r c 0o)]- Let us 
point out also the high energy regime corresponds to a 
4d regime; i.e. H 2 ~ n\p/Z. To complete our analysis, 
based on Eqs. ([3TT]) . (pH|) . (j3~3"| and (JSU), is it is helpful 
to distinguish two cases: a positive and a negative 70. 



1. Positive 70 

The solution Z + corresponds to a brane that starts 
its evolution with a Big Bang singularity where p ~ 
4/(37oK|)i -2 . Consequently, at early time, coinciding 
with the high energy regime, the scale factor scales as 
a ex i 2 /( 3 T°) corresponding to the relativistic behaviour 
of a 4d Friedmann-Lemaitre-Robertson- Walker (FLRW) 
universe filled with the same matter content as the brane. 
This is in agreement with the fact that a homogeneous 
and isotropic DGP brane at high energy behaves like a 
FLRW universe in 4d general relativity pj], EI- The 
brane keeps expanding and at very late-time, the en- 
ergy density and scale factor can be approximated by 
p ~ l/(ft!r c 7o)£ -1 and a ex t 1 ^ 3 ' t0 \ respectively. This is 
a consequence of the low energy behaviour of the nor- 
mal DGP branch (with a + sign in Eq. (12.31) ) where 
H~r c K\p/2, Pill. 



The solution Z_ is asymptotically de Sitter on the past 
where p ~ exp (37o£/ r e) and a ex exp(— t/r c ). Notice that 
t is negative as can be easily seen from Eq. (|3.2p bear- 
ing in mind that 70 is positive. Afterwards the brane 
starts contracting until it hits a big crunch in the future, 
where p ~ 4/(37^|)t _2 and a ~ t 2 /07o)_ Because the 
brane is contracting in this case and < 70, the high en- 
ergy regime takes place at late-time while the low energy 
regime describes the early-time evolution of the brane. 



2. Negative 70 

For a negative 70, the energy density is a growing 
function of the scale factor. Therefore, the high-energy 
regime corresponds to large scale factors while the low- 
energy regime takes place at small scale factors. 

The brane described by Z + starts its expansion at 
t — > — 00 where a ~ 0. More precisely, at early time 
p ~ l/(K 2 r c 7o)£ _1 and a 370 ex jot. This regime cor- 
responds to the low-energy regime of the normal DGP 
branch (with a + sign in Eq. (|2.3[1 ) where the Hubble rate 
is approximately a linear function of the energy density 
of the brane. The brane keeps expanding in a super- 
accelerating way; i.e. < H, as Eq. (|2.6p implies for 
70 < 0, < p and e = 1. The expansion of the brane halts 
when it hits a big rip singularity at t = where the scale 
factor, the energy density, the pressure, the Hubble rate 
and its cosmic derivative blow up [22| . |23[ . The last stage 
of the brane expansion corresponds to the high energy 
regime; therefore a 4d regime, where p ~ 4/(37qk|)£ -2 
and a 37 ex t 2 while t — > 0~. 

On the other hand, the brane corresponding to Z~ 
starts its evolution with a type I singularity [23| in the 
past. Initially, the energy density and scale factor are 
very large, indeed p ~ 4/(37gK 2 )£ -2 and a 37 ex t 2 
where t — > + . This corresponds to the high energy 
regime. Then the brane starts contracting; i.e. H < 0, 
even though < a because < H. The brane is 
asymptotically dc Sitter in the future, corresponding to 
the low energy regime, where p ~ p m exp(37oi/ r c) an d 
a ex exp(— t/r c ). 



B. Negative energy density 

The DGP model has solutions with a A:-field with a 
finite negative energy density as long as — p m < p < ; 
i.e. Z < -1, cf. Eqs. (H3J) and ([2TT2|) . At the mini- 
mum energy density p — — p m , which corresponds to a 
finite value of the scale factor, the Hubble rate vanishes 
(see Eq. (|2.3p ). while the derivative of the Hubble rate 
diverges (cf. Eq. (|2.6p ). eventhough the energy density 
and the pressure are finite. Therefore this is a "quies- 
cent" singularity [241 . but notice that it is different from 
the one discussed in [24| which was induced by: 

• the presence of a dark radiation term on the brane, 
and therefore a black hole in the bulk which is not 



the case in the DGP model, and matter with pos- 
itive energy density and whose equation of state 
satisfies specific conditions. 

• or through an inequality condition satisfied by a lin- 
ear combination of the bulk cosmological constant 
and the brane tension which cannot be fulfilled on 
the DGP model. 

Note as well that a "quiescent" singularity can also shows 
up in a DGP brane-world model with a Gauss-Bonnet 
term in the bulk [25|. Here again, the singularity is dif- 
ferent from the one that takes place at p — —p m . The 
singularity in [251 ] is caused by a combination of infra-red 
and ultra-violet modifications of general relativity. The 
presence of this singularity is the main motivation for 
looking for the solution of Eq. (|2.12[) when the energy 
density is negative. We summarise next the dynamics of 
the normal branch on those cases. 

For p < 0; i.e. Z^ < 0, Eq. (|2.12[) can be rewritten as 
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(3.5) 



which can be also integrated analytically [2C 
1 
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(t-t 2 ), 



(3.6) 



where t 2 and K 2 are constants. Here again, the function 
Z + , Z~ will refer to the Z function satisfying Eq. (|2.12[) 
with + and — , respectively. 

For p < with p < 0; i.e. Z ± < 0, the solution (ESJ) 
can be expanded about the minimum energy density, p = 
-p rn or Z = -1, 



37o i± 2 
2r n 3 



37o , 
2r c 



(3.7) 



where a^J° = —po/p m - Below, we will see that the sec- 
ond term in the expansion is responsible for the diver- 
gence of a when it is evaluated at a = a m or t = 0. 
At this point the brane approaches the "quiescent" sin- 
gularity because H remains finite but H diverges. For 
simplicity, we have set a = a rn at t = 0. The po- 
tential close to the singularity can be approximated by 
V ± ~ V a-?~< [1 - 37o0/(2r c o )]- 1 and it is clearly fi- 
nite at the singularity where = 0. The scale factor a m 
stands for the scale factor at the "quiescent" singularity 
which is the minimum scale factor of the brane for < 70 
and the maximum scale factor of the brane for 70 < 0. 

On the other hand, at very low energy (p — > 0~ or 
Z — > —00), we obtain 



Z+ 

H-z-) 



3to j 
4tv 
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(3.8) 
(3.9) 



The /c-field potential on this regime fulfils: V + ex 4>a/cj) 
and V~ ex exp[37o0/(r c 0o)]- It is also worthy to mention 
that the brane always shrinks when filled by a negative 
energy density (see Eq. (|2.3p with e = 1). Like in the 
previous subsection, to complete our analysis, based on 
Eqs. (|3TTj) . (j3~7| . (|3~5|) and ([SI]), is it is helpful to distin- 
guish two cases: a positive and a negative 70. 



1. Positive 70 

For Z + and < 70, the brane starts its evolution with 
a very large radius and an almost vanishing energy den- 
sity, where p ~ por c /t and a ex (t/Vc) 1 / 370 . The Hubble 
rate and the derivative of the Hubble rate are almost 
zero initially. This solution is defined in a region where 
the cosmic time is negative. Afterwards the brane starts 
contracting until it hits the "quiescent" singularity at a 
minimum scale factor 



37 t 2 / 37 t 



a ~ a r , 



2r r 



t 2 

2V C ~ 9^ 
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, (3.10) 
, (3.11) 



where the latter equation shows explicitly that a diverges 
at t = 0. 

The solution corresponding to Z~ is asymptotically de 
Sitter in the past where once more the scale factor is 
very large and the energy density is approaching zero. In 
more specific terms, we obtain p ~ — p m cxp(37o£/r c ) and 
a ex exp(— t/r c ) where t <C — 1. Even though the brane 
is asymptotically de Sitter in the past, it starts shrink- 
ing because the Hubble rate is initially (almost) constant 
and negative. In the future the brane face a "quiescent" 
singularity at a minimum scale factor characterised by 



37 t 2 / 37 t 
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, (3.12) 
, (3.13) 



and again a diverges at t = 0. 



2. Negative 70 



The solution Z^ starts at a "quiescent" singular- 
ity, where the brane reaches its maximum scale factor 
a m , and has the same behaviour described by the so- 
lutions (|3.10p - (|3.13j) . In this case the solutions are de- 
fined in the region where the brane cosmic time is pos- 
itive. The solutions Z^ describe a contracting brane 
that ends with a vanishing scale factor even though at 
different rates. For Z + , the energy density and scale 
factor change with the cosmic time as p ~ por c /t and 



a ex (t/rc) 1 / 370 , while for Z~, the same quantities be- 
haves as p ~ — / o m exp((37o£)/r c ) and a oc exp(— t/r c ). 
In both cases, the brane shrinks to a point in an infinite 
amount of its proper time. 



IV. CONCLUSIONS 

We have investigated a DGP brane filled with a k- 
cssencc field. This source is particularly interesting be- 
cause it can take negative values and gives the possibil- 
ity of considering "quiescent" singularities in the DGP 
framework. In addition, it gives rise to inflationary 
branes. We have focused our analysis on the normal 
DGP branch, for the reasons explained in the introduc- 
tion, fundamentally based on the duality between both 
branches as we have shown at the end of section II . Our 
results for the normal branch are summarised on the next 
paragraphs. 

In the < p case we have found: (i) a singular scale 
factor that represents an expanding universe with a final 
power law scenario and, a contracting one with an initial 
de Sitter phase and a final big crunch at t = for 70 > 0. 
However, for 70 < 0, (ii) we have found an expanding 
scale factor, such that, the universe begins with a van- 
ishing value at t = — 00 and ends in a big rip singularity 
at t = 0. There exists also a singular universe which 
contracts from a = 00, at t = 0, to an asymptotically de 
Sitter scenario in the future. 

In the p < case we have obtained a contracting scale 



factor for any value of the equation of state 70. For 70 > 
the universe contracts from the far past at t = —00 with 
a power law expansion or a de Sitter phase and ends with 
a constant scale factor a m , at t = 0. We have shown that 
the expansion rate is finite; more precisely H = H m = 
— l/2r c , at t = and H m diverges, indicating that the 
universe ends in a quiescent singularity [241 ] . Nevertheless 
the potential becomes finite when it is evaluated at this 
singularity, V m = Vaa^ 10 . Similarly, the energy density 
and the pressure of the /c-field are finite at this point. For 
7o < the universe begins from the constant value a m 
in a "quiescent" singularity, after that, it contracts and 
ends with a power law or a de Sitter phases. 
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